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Abstract 



The results of ab initio calculation of El amplitudes, lifetimes, and polariz- 
abilities for several low-lying levels of ytterbium are reported. The effective 
Hamiltonian for the valence electrons H e q has been constructed in the frame of 
CI+MBPT method and solutions of many electron equation Pl c ff*& n — E n *§> n 
are found. 
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I. INTRODUCTION 



In this paper we report results of an ab initio calculation of El amplitudes, lifetimes, and 
polarizabilities for several low-lying levels of ytterbium. In Ref. [Q we calculated the energies 
and hyperfme structure (hfs) constants of low-lying levels of ytterbium. In that calculations 
the accuracy of atomic wave functions was tested at the short distances by comparison of 
the calculated hfs constants with the experimental ones. The latter are usually known to 
a very good accuracy, providing a good test of the quality of the wave function near the 
nucleus. 

El amplitudes, in contrast, are determined by the wave function behavior at large dis- 
tances. Usual experimental accuracy for the oscillator strengths and scalar polarizabilities 
is on the level of few percent. This is close or even less than the accuracy of precise atomic 
calculations (see, e. g., calculations for Ba and Cs ||). Tensor polarizabilities can be 
measured with the accuracy of 1% or better Thus, it is possible to test an atomic wave 
function at large distances at 1% level. Note that 1% accuracy is crucial for calculations of 
parity nonconservation effects in atoms, because it allows to test predictions of the Standard 
model at small momentum transfer ||||. So far such precision has been achieved only for 
one-electron atoms Cs and Fr |7|-|9| • I n this work we deal with a much more complicated Yb 
atom. 

We consider ytterbium as a two electron atom with the core [Is 2 , 4/ 14 ]. Valence- 
valence correlations are taken into account by the configuration interaction (CI) method, 
while core- valence and core-core correlations are treated within the second order of the many- 
body perturbation theory (MBPT). The latter is used to construct an effective Hamiltonian 
for the CI problem in the valence space. The details of the method can be found in the 
papers [|1^,|11|]. Application of this method to calculations of hfs constants has been discussed 
in |Jiy2],|12"|l . In Ref. [@ the method has been extended to the calculations of polarizabilities. 
Here we apply this technique for calculations of lifetimes, El amplitudes, and polarizabilities 
of ytterbium. 

II. GENERAL FORMALISM 

Let us write out several formulae that will be used in the following. The expression for 
oscillator strength for a, J —>■ a', J' transition has the form |13| (atomic units m — % — e — 1 
are used throughout the paper): 

f{ala'J') = -^f Y) \{a,J\\D\\a\J')\\ (1) 

where {jJ a j,a<j< = E aJ — E a >ji, D is the dipole moment operator, and reduced matrix elements 
(MEs) are defined as follows: 

(a',J',M'\D q \a,J,M) 

= (-!)'-"' ( _ J M , \ J M )«f\\D\\a,J). (2) 

The lifetime r of a level is the inverse of the total transition rate. The probability for 
a, J — > a 1 , J' transition is given by: 
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W(aJ, a' J') = ^f^\{a, J\\D\\a>, J')\\ (3) 

where c is the speed of light. 

Static polarizability of the sublevel \a, J, M) in a DC electric field E = £z is defined as: 

A£ M , A/ = --a a ,j jM £ 2 (4) 

if 3M 2 -J(J+l)\ c2 
- — 2 [a , a ,j + «W J{2J _ l) ) £ > 

where AE a j t M is the energy shift and qjq and ct2 define the scalar and tensor polarizabilities, 
correspondingly. Being a second order property, a a ^ t M can be expressed as a sum over 
unperturbed intermediate states: 

_ | (a, J,M\D z \n, J n ,M)\ 2 
a a ,j,M = _ , (5) 

where E n is an unperturbed energy of a level n, and the sum runs over all states of opposite 
parity. The formalism of the reduced MEs allows to write explicit expressions for the scalar 
and tensor parts of the polarizability: 

-2 \(a,J\\D\\n,J n }\ 2 

«o,M- 3(2J+1) 2 r K _ En , W 
/ 40J(2J-1) V /2 



i 3(2J + 3)(2J+l)(J+l) 

^ J+J „ +1 / J 1 J n \ \(a, J\\D\\n, J n )\< 
1 J 2 — E„ 



E(-i) 



In order to use Eqs. ©-(0) in calculations one needs to know a complete set of eigenstates 
of the unperturbed Hamiltonian. It becomes practically impossible when dimension of a CI 
space exceeds few thousand determinants. It is known, that it is much more convenient to 
solve inhomogeneous equation instead of the direct summation over the intermediate states 



14]|L5|1 - Indeed, let us consider the solution of the following equation: 

(E a -H)\X aM ,) = D q \a,J,M), (8) 

where q = 0, ±1 and M' = M + q. Obviously, the right hand side in Eq. (|5|) can be expressed 
in terms of the function X ai M (note that D Q = D z ): 

a a ,j,M = -2 (a, J, M\D \X aM ). (9) 

If we want to rewrite Eqs. @ and (^) in terms of the function X Qj m', we need to de- 
compose the latter in terms that correspond to particular angular momenta Jj. Generally 
speaking, there can be three such terms with Jj = J, J ± 1: 

X a ,M' = -Xo.J-l.Af' + X a> j t M' + X a> j + i t M'- (10) 
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Now, with the help of the functions X^j^m' Eqs. (JD) and (^) are reduced to: 

a , a ,j = (-l) g+1 Q , OT 2 , r , (11) 

2 

(a, J,M\D_ q \X a ,j>, M ,), 




V. 




1/2 




3)(2J + 1)(J 

j i j' r 2 

1 J 2 J 

-2 

(a, J,M| J D_ 9 |X a ,j/ i Af), 



(12) 



where sums run over J' = J, J± 1. Note, that these equations are valid only if all 3j-symbols 
on the right hand side do not turn to zero. One has to take it into account when choosing 
for what spherical component q to solve Eq. (|8|). 

If we know the solution of Eq. (|j) and its decomposition ([U]), then expressions (|11|) and 



(|12|) allow us to find both scalar and tensor polarizabilities of the state \a, J). Moreover, the 
same functions X a> ji y M' can be also used to find other second order atomic properties, such 
as amplitudes of the Stark-induced El transitions or parity nonconserving El transitions 
between the states of the same nominal parity (see, for example, Ref. ||16||). 

III. CALCULATION DETAILS AND RESULTS 

A. Orbital basis set and CI space. 

The calculation procedure is quite similar to that of Ref. 0. For this reason we give here 
only a brief description of its features. This calculation is done in the V N approximation, 
that means that core orbitals are obtained from the Dirac-Hartree-Fock (DHF) equations for 
a neutral atom (we use the DHF computer code [ffTJ). The basis set for the valence electrons 
includes 6s, 6p, 5d, 7s, 7p, 6d DHF orbitals and 8s-15s, 8p-15p, 7d-14d, 5f— lOf, and 5g-7g 
virtual orbitals. The latter were formed in two steps. On the first step we construct orbitals 



with the help of a recurrent procedure, which is similar to that suggested in Ref. |L8[ and 
described in Refs. [pl|IB|. After that we diagonalize the V N DHF operator to obtain the 
final set of orbitals. 

For this orbital basis set the complete CI is made for both even-parity and odd-parity 
levels. Two-electron wave functions are the linear combinations of the Slater determinants 
with a given J z . It means that no symmetrization with respect to angular momentum J is 
made. 



B. Effective operators. 

Within the CI+MBPT method the wave function of the valence electrons is found from 
the eigenvalue equation: 
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H eS \a, J,M) =E a \a, J,M). (13) 

Eq. (H) is rewritten as equation for valence electrons only: 

(E a - H eS )\X a>M ,) = D cff >, J, M>, (14) 

with the effective operators, which are found by means of the MBPT. The effective Hamilto- 
nian for two valence electrons is formed within the second order MBPT []TD1 . We used RPA 
for the effective dipole moment operator (see, for example, Ref. |[19|| ). We have checked that 
MBPT corrections to D e g, which are not included in RPA, are small if RPA equations are 
solved with 6s electrons excluded from the self-consistency procedure. That means that RPA 
equations have the same form as in the V N ~ 2 approximation. The more detailed description 



of the effective operator formalism is given in [12 



C. Transition amplitudes and lifetimes. 



We first solve eigenvalue Eq. fll3|) with the effective Hamiltonian for low-lying even- and 
odd-parity states. Strictly speaking, the effective Hamiltonian can be safely used only for the 
energy levels below the core excitation threshold. For Yb this threshold lies at 23189 cm -1 



above the ground state |20]. However, it was shown in [|TJ that theoretical spectrum is 
quite good up to ~40000 cm -1 . Correspondingly, we can work (with some caution) with 
the states lying slightly above the core excitation threshold. In our approach we fail to 
reproduce the states with unfilled / shell and correspondingly to account properly for the 
interaction with such states. For this reason we restrict ourselves to the consideration of 
the states lying sufficiently far from those with unfilled / shell. We consider El transitions 
between four low-lying odd-parity states ( 3 Pq 12 (6s6p) and 1 P° (6s6p)) and seven even- 
parity states ( 1 S'o (6s 2 ), 3 -Di,2,3 (5d6s), X D 2 (5d6s), 3 Si (6s7s), and 1 So (6s7s)). The state 
1 P° (6s6£>) requires special attention. The nearest J 13 5d 6s 2 state lies only 3800 cm -1 above 
the latter and their interaction is not negligible. We estimated that configuration / 13 5<i6s 2 
contributes on the level of several percent to the wave function of l P° (6s6p) state. We do 
not take into account this configuration mixture. This reduces the accuracy of the calculated 
1 P° (6s6p) — > 1 Lj El amplitudes. 

When eigenfunctions for the valence electrons are found, we can calculate transition am- 
plitudes and lifetimes. The results of calculations are presented in Table |. The magnitudes 
of the El amplitudes vary in a wide range. These variations correspond in part to the 
approximate selection rules AS = and A J = AL, which are easily traced through Table |. 
For large amplitudes we estimate the accuracy of our calculation to be 3-5%. For the reason 
discussed above the amplitudes ( 1 Lj|D| 1 P 1 ° (6s6p)) do not follow this rule. The accuracy 
for these amplitudes, as well as for small amplitudes (< 0.5 a.u.), is about 15-20%. 



Where available, we compare our results with those of other theoretical [plp2[ and ex- 
perimental P3"-28| groups. For the convenience of comparison we recalculated the oscillator 
strengths and transition probabilities to the reduced MEs. Calculations in Ref. |21j were 



performed in the L-S coupling scheme. The simplest semiempirical method p9| was used 



then to evaluate the radial parts. In Ref. [22] the mufticonfiguration Dirac-Fock method 



was used. The valence-core electronic correlations were included semiempirically. Compar- 
ing our results with the results of other theoretical works, one can see that it was important 
to account for the valence-core correlations. 
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Now, using Eq. (|3]) we can find the transition probabilities and the lifetimes of the levels 
(see Table |T|). In these calculations we used experimental transition frequencies. Therefore, 
the accuracy of these numbers depends only on the accuracy of the dominant transition 
amplitudes. As a result, the largest error (40%) takes place for the states 1 S , (6s7s) and 
1 D2(5d6s) where the transition to the state 1 P 1 °(6s6p) is dominant. For other states we 
estimate theoretical accuracy for the lifetimes as 10% or better. 



D. Polarizabilities. 

In order to find the polarizabilities we substitute eigenfunctions in the right hand side 
of Eq. (0) and solve corresponding inhomogeneous equation. After that Eqs. (|TT1) and ( |12|) 
give us a and a 2 . Results of these calculations are presented in Table [TTT|. It is seen that, 
unlike barium (see Ref. @), a 2 has typically the same order of magnitude as a . For this 
reason the theoretical accuracy for a 2 , as a rule, is similar to that for a - I n contrast, 
experimental data for a 2 are usually much more precise and complete. 

There are several sources of errors in the calculations of polarizabilities. Some of them 
are the same as for hfs calculations, and are connected with the inaccuracy in the wave 
functions and the effective operators (note, that RPA corrections to the dipole operator are 
much smaller than for hfs operators). The additional source of errors is the inaccuracy in 
eigenvalues. Finally, solving Eq. (|i~4l) we do not account for configurations 4/ 13 nl n'V n"l" . 
In other words, we do not take into account / shell polarization. Fortunately, the states 
of configuration 4/ 13 nl n'V n"l" , that can be reached by one-electron transitions from the 
levels studied here, lie rather far above the latter. The estimates show, that contribution of 



/ shell polarization to the polarizabilities of the states listed in Table [III] does not exceed 
2-3 a.u. 

The final accuracy of calculations is very different for different levels. For instance, the 
95% of the polarizability of the ground state 1 S (6s 2 ) is due to the ME ^SqID^P" (6s6p)). 
Supposing that this ME is calculated within the accuracy of 20%, the latter for a ( 1 5 , ) 
will be about 40% (corresponding transition frequency is reproduced almost ideally |J). We 
need to say that even taking into account the large uncertainty of our result, it significantly 
differs from «o ( 1 *5'o) = 266 a.u. obtained in Ref. where the Hartree-Fock method was 
used. 

For the Dj (5d6s) states the situation is more complicated. There are large cancellations 
between contributions of Pj (6s6p) states and higher-lying states. For this reason their 
polarizabilities are small and the role of different small contributions is enhanced. Thus, 
analysis of the accuracy becomes difficult; only for the tensor polarizability of 3 D^ (5d6s) 
state we can estimate the accuracy to be 20%. All other values of «o and a 2 for Dj (5d6s) 
states presented in Table fTl] are rather estimates by the order of magnitude. 

The scalar polarizabilities of the levels 3 Si (6s7s) and 1 S , (6s7s) are basically determined 
by the MEs (Sji (6s7s)\D\Pj (6s7p)). Because of the closeness of / 13 5d 2 6s states we failed 
to obtain the reliable wave functions for Pj (6s7p) states. Correspondingly, the values for 
«o( 3 <Si (6s7s)) and «o( 1 5'o (6s7s)) are also only the estimates. 

Now let us go over to the odd-parity states. The accuracy of «o and a 2 for 3 Pj (6s6p)- 
triplet is 6-10%. The main contribution here comes from the 3 Dj (5d6s) multiplet and there 
are no cancellations because all important levels of opposite parity lie above and contribute 
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with the same sign. The accuracy for a of 1 P°(6s6p) state is about 40% and for a 2 even 
worse (50%). This is due to the large contribution of the intermediate state 1 So (6s7s) to 
these polarizabilities (see above). 

In Ref. J3l| the Stark shift of the 1 5'o (6s 2 ) — > 3 P° (6s6p) transition in ytterbium was 
measured. The Stark shift rate was found to be K= —61.924(0.193) a.u. In terms of 
polarizabilities this magnitude can be written as: 



K = -\ {«o ( 3 A°) " 2a 2 ( 3 iT) - a CS )} ■ 



Using the numbers from Table |T|, we find that K= —55 (9) a.u., in good agreement with 



the experimental result pi 



In Ref. ||23|| , the Stark shifts for x Sq (6s 2 ) — > 3 Di^ (5d6s) transitions were observed. These 
shifts depend on the differences in scalar polarizabilities (ao ( 1 5'o) — «o ( 3 -Di,2)): 

J 71 theory, 
1 86 (3) experiment, 



a CS )-a CD 2 ) = lll (A , the ° r ^ 
u v u; u v 21 [80 (4) expeni 



experiment, 

where theoretical values are taken from Table f[\\ . 

The method used here allows us to calculate not only static polarizabilities, but, also 
the Stark-induced amplitudes for different transitions. For instance, the magnitude of the 
vector transition polarizability \/3\ for 1 S'o (6s 2 ) — ► 3 Di (5d6s) transition was calculated to 



be 122(12) a.u., in good agreement with our previous calculation 138(30) a.u. p2| and 



experimental result 113(15) a.u. [23] 



E. Conclusion. 

Application of the effective operator technique to the Yb atom is hampered by the ex- 
istence of the shallow 4/ shell. Nevertheless, it is possible to make reliable calculations of 
different atomic properties including transition frequencies, hyperfine constants, El ampli- 
tudes, lifetimes and polarizabilities for many low-lying energy levels. It is of a particular 
importance, that with some caution calculations can be done even for levels above the core 
excitation threshold, which for Yb lies at 23189 cm -1 . 
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TABLES 



TABLE I. Reduced MEs |(-Lj||r||Xj,)| (a.u.). Calculations were made in the L-gauge. Other 
theoretical and experimental results are given where available. The uncertainties are indicated in 
the parentheses. 





d P °(6s6p) 


^P°(6s6p) 


A Pl ?(6s6p) 


L P°(6s6p) 


l S {6s 2 ) 





0.54(8) 





4.40 (80) 






0.44 a 




4.44 a 






0.549 (4) c 




4.89 b 






0.553 (13) d 




4.13 (10) d 










4.02 c 










4.26 f 


J-/\ \Oil\Jo j 


Z.Ul y 1U J 


9 9fi Hfll 


U.UU ^ X.Z ) 


n 97 (~\ (\\ 






2.2(1) § 




0.24 a 


3 D 2 {5d6s) 




4.03 (16) 


2.39 (10) 


0.32 (6) 










0.60 a 


3 D 3 (5d6s) 






6.12(30) 




1 D 2 (5d6s) 




0.54(10) 


0.38 (8) 


3.60 (70) 


3 5i(6s7s) 


1.98(10) 


3.53(15) 


5.05 (20) 


0.73(15) 




1.36 b 


2.50 b 


3.77 b 




1 S* (6s7s) 




0.22 (4) 




4.31 (80) 






0.22 (2) h 






Theory: a Ref. [p2|, 
Experiment: c Ref. 


b Ref. 0; 
23 , d Ref. 


24 , c Ref. 25], f Ref. 26], 


g Ref. [27 , h Ref. [28]. 





TABLE II. Lifetimes (nsec) of the low-lying 


levels for Yb calculated with the reduced MEs 


from Table [ 


] and experimental transition frequencies. 




State 


Config. 


This work 


Other data 


3 D 1 


5dQs 


372 (30) 


380 (30) a 


3 D 2 


5dQs 


430 (35) 


460 (30) a 


3 D 3 


5dQs 


540 (55) 




l D 2 


5dQs 


4400 (1800) 


6700 (500) a 


3 s 1 


6s7s 


13.5(1.1) 


12.5(1.5) b 








15.9(1.9) c 


'So 


6s7s 


33(13) 


45.8(1.0) d 


3po 


6s6p 


875 (250) 


760-875 e 








1294 f 


3 no 
r 2 


6s6p 


15.0(1.5) sec 


14.5 sec f 


lpo 


6s6p 


5(2) 


5.1-6.4 e 



4.78 



Experiment: a Ref. ]^7||, Ref. []3"3" |, c Ref. [S3], Ref. |[28|| , c see Ref. [3_S] and references therein; 
Theory: f Ref. [||]. 
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TABLE III. Scalar and tensor polarizabilities (a.u.) of low-lying levels of Yb. Theoretical 
accuracy is indicated where analysis was possible, otherwise the numbers should be considered as 
estimates. 



Theory Experiment 



level 


config. a 


«2 


«2 


'So 


6s 2 118(45) 






3 D 1 


5d6s 47 


22 


28 (4) a 


3 D 2 


5d6s 36 


17 


28 (8) a 


3 D 3 


5d6s -9 


118(24) 




X D 2 


5d6s 4 


150 






fi„7„ 9D79 
OS 1 S ZU 1 z 








6s7s 2030 


0.8 




3 no 
-M) 


6s6p 252 (25) 






3po 


6s6p 278 (15) 


24.3(1.5) 


24.06(1.37) b 
24.26 (0.84) c 
23.33 (0.52) d 


3 DO 

r 2 


6s6p 383 (32) 


-76 (6) 




lpo 


6s6p 501 (200) 


-118(60) 


-57.4(5.6) b 


a Ref. [23], 


b Ref. H, c Ref. 37 , d Ref. gl] . 
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